On the conversion of mass eigenstates 



Mikhail V. Medvede\H 

Institute for Advanced Study, Princeton, NJ 0854 

In this paper we consider a stable particle with flavor mixing. We demonstrate that incoherent 
conversion of heavy mass eigenstates into light ones and vice versa can occur, as a result of elastic 
scattering. This effect is nontrivial for non-relativistic particles, for which the standard flavor 
oscillation ceases rapidly due to incoherence. We also prove that if a heavy state is bound in a 
gravitational potential and a light state is unbound, the mass-state conversion can lead to gradual 
"evaporation" of the mixed particle from the potential. A number of implications, ranging from the 
cosmic neutrino background distortions to scenarios of cold dark matter evaporation from halos, are 
addressed. 



PACS numbers: 03.65.-w, 14.60.Pq, 14.80.-j, 95.35.+d 

Introduction — Massive flavor-mixed particles exist: 
the examples are quarks, K°s, neutrinos, etc. Dark mat- 
ter can also belong to this class. The propagation (mass) 
and interaction (flavor) eigenstates of mixed particles [1 
are related by a unitary transformation, 

l/i) (1) 

3 

where |/) and \m) denote the flavor and mass states, 
and U is a unitary matrix. If the above superposition is 
coherent, then the time-dependent interference of mass 
states propagating with different speeds leads to the fla- 
vor oscillation phenomenon, which has been experimen- 
tally observed for relativistic neutrinos, for instance. In 
contrast to flavor eigenstates, mass eigenstates are more 
fundamental: they are thought to exhibit no transfor- 
mations (except for quantum broadening of their wave 
packets and/or decay into lighter particles if the primary 
is unstable) and are usually treated as "real particles" in 
a semi-classical sense. Indeed, the mass state wave pack- 
ets have, in general, different energies and momenta and 
they propagate in the same way as the normal particles 
of masses would I2j , whereas the observed pure flavor 
states are merely a time-dependent interference pattern. 
This view is also supported by the fact that it is impos- 
sible to construct a Fock space of flavor states, yet flavor 
oscillations can be derived without the concept of flavor 
eigenstates [3]. 

In this paper we show that transformations (or conver- 
sion) of mass states, \irii) — > J2j a j \ m j) ( a j being com- 
plex amplitudes), are possible even for stable particles. 
Repetitive conversions, schematically shown as 

... -> m h -> m h -> m h -> . . . 

\ \ \ (2) 
mi mi 

can have interesting implications for astrophysics and 
cosmology. We will make a number of assumptions 
throughout, which are useful for presentation purposes 
but are not crucial for the phenomenon at hand, hence 
they can be generalized or omitted. 



Description — Let us consider a two-flavor system 
of stable particles, for simplicity. Without loss of gen- 
erality, mi > m,%. If m\ ^ m%, we refer to them 
as the heavy and light states, \mh) and \mi), respec- 
tively. The interaction matrix, diagonal in the flavor ba- 
sis, V — diag(V r a , Vp), is non-diagonal in the mass basis, 
V = IPVU. In particular, the off-diagonal terms are 
Vlh = Vm = (V a — Vg)cos#sin#, where is the mixing 
angle and we assume V a ^ Vp. 

It is convenient to use the wave packet treatment, 
in which mass states are gaussians whose peaks follow 
classical trajectories. In general, mass states propagate 
with different velocities and can separate in the spa- 
tial domain. Such a mass state separation means that 
the wave packet overlap diminishes and the particle's 
wave-function becomes de-localized (not singly-peaked in 
space). Hence, the superposition of mass states is no 
longer coherent, so that the probability of detection of a 
certain flavor ceases to oscillate and tends to a constant, 
set by the mixing matrix. Along with the described bal- 
listic mass-state separation, they can also separate via 
the (mass-dependent) interaction with gravitons or via 
their motion along different geodesies in a static gravity 
field (see more discussion below). 

As an illustrative model, let us consider a particle in a 
one-dimensional box. The left and right ends of the box 
are "membranes" (or semi-transparent "mirrors") that 
keep the heavy states inside but let the light states to 
freely escape from the box, see Fig. [l] Somewhere in the 
middle of the box, we put an elastic scatterer with an 
interaction matrix, V. We assume that the creation, de- 
tection and scattering processes are such that the widths 
of the mass state wave packets are always smaller than 
the box dimensions; this can be done by choosing a suffi- 
ciently large box. Now, consider a particle that bounces 
back an forth in the box and experiences elastic scat- 
tering from time to time. For generic initial conditions, 
the membranes let the light states escape from the box, 
so upon one transit time the box will contain a heavy 
state only. Because of off-diagonal terms in V, scatter- 
ing of a single mass state yields a mixture of both states 



2 



1 

M 






M 






J a, \m,)+a h \m h ) 


a, |m,>^ 


1 
1 

a',\m,) | 


f 5 * " 

L_ a',|m ; >+a',,|m,,> | 


I |m,,) 






1 
1 









FIG. 1. Cartoon of a mixed particle in a box. Left and right 
ends of the box are "membranes", M, which keep a heavy 
state inside and let a light state to escape. A block in the 
center depicts a scattering potential V . 

with smaller amplitudes. (Note that scattering shall not 
happen too frequently, because the mass states have to 
be well separated; otherwise the mass states may inter- 
act coherently, as a well defined flavor state.) Further 
evolution and interaction with a membrane leaves only a 
heavy state in the box again. Thus, by the end of the 
cycle, the system contains the same mass state but with 
a smaller amplitude. By repeating the cycle many times, 
one effectively converts the heavy mass state in the box 
into the light one represented by a "wave-packet train" 
outside the box, as in Eq. (|]). In a sense, the parti- 
cle "evaporates" from the box. We are not aware of any 
discussion of such a process in the literature. 

This process can be described better using the oper- 
ator formalism. The processes of elastic scattering and 
mass state separation due to propagation and interaction 
with a membrane can be formally described by operators 
V and Ai, which transform the mass state amplitudes 
and coordinates, respectively. Since the scattering po- 
tential is localized within the box, it has a compact sup- 
port, Supp(F) C Lbox, where Lb ox = {1 | < i < 
ibox} C K is a one-dimensional compact manifold rep- 
resenting space inside the box of length Lbox- There- 
fore, we define the V-operator as V: {a,i} —> {o^} such 
that V J2i a t( x i)\ m i) = Ei a 'i( x i) \ m i) > i = l,h with 
a'^Xi) = VijCLj(xj), Xi = x.j Mxj e L box and a'^xi) = 
ai{xi) yxi (jL Lbox, where the transition matrix elements 
Vij oc Vij computed on the shell ("hat" denotes tran- 
sition amplitudes to distinguish them from the matrix 
elements of the scattering potential), a,i(xi) are complex- 
valued functions (representing shape factors and ampli- 
tudes of wave packets) of one argument Xi , respectively. 
The VW-operator is defined as M. : {ai\ — > {a.i} such 
that M J2i a i( x i)\ m i) = J2i a i( x 'i) l m ») . xi,x h ,x' h € 
L b ox, x'i G ^\Lbox- Generalization to three dimensions 
and more than two flavors is straightforward. For the 
sake of simplicity, we also assume that the wave pack- 
ets outside the box remain localized and never overlap 
with each other, as well as never come back into the box. 
We also simplify Vij by assuming that one of the fla- 
vor states is non-interacting, Vg = 0, and another one is 
strongly interacting (with unit probability), V a = 1, and 



non-vanishing mixing is also assumed. Generalization of 
these assumptions is straightforward as well. 

The following theorem holds true: Let 'Ffc = 
EiafVi). k = 0,1, then Va, (0) : Supp(af ) ) C 

Lbox, E>f| 2 = hj = l>h> 3a« : EJ^I 2 = 1 
that satisfy = XiVAi 'Fo, and the following is true: 

l4 X) | 2 < l4 0) | 2 . K (1) | 2 > K (0) | 2 ^d Supp(4 1) ) C Lbox, 
where the norm is defined in a standard way: |ai| 2 = 
/ a*a,idxi. That is, for an arbitrary initial state, <F, of 
a mixed particle which is initially in the box, the total 
detection probability of the heavy state decreases upon 
the action of the composite operator MVM. upon 'F, yet 
it remains identically zero outside the box. 

The proof follows from the direct substitution. First, 
let's label coordinates inside and outside the box as x\ £ 
Lbox and x[ € K\Lbox- Next, we have 

*i = MVM 2 4°\xj)\mi) 

i—l,h 

= MV (af\x\) \mi) + a^(xi) \m h )) 

= M (af\x\) \mi) + af\xj) |m,) + a'^(xi) \m h )) 

= (a^ixf) + af\xj)) |m,) + af\xi) \m h ) . (3) 

Thus, = a^'Vhh and = +a^Vih- Obviously, 
[a^j 2 = |a^| 2 |Vftfe| 2 < \ a ^\ 2 f° r any nonzero mix- 
ing. Also, h (1) | 2 = |a| 0) | 2 + |ai 0) | 2 |Vk| 2 > |a| 0) | 2 , where 
we used that wave packets do not overlap. Obviously, 

E, l^l 2 = l«| 0) | 2 + l4 0) | 2 (l%| 2 + \%h?) = 1. Finally, 
from eq. (J3J) , one has Supp(a| 1 ' ) ) = Supp(a'[°' ) ) C Lbox- 
q.e.d. 

Clearly, multiple application of the (A'fV)-operator re- 
sults in further reduction of the heavy state probability 
amplitude. If we define (MV) n M = (MV) . . . (MV)M, 
then it is straightforward to show that for >F„ = 
(MV) n M$o one has \a^ ] \ 2 = \ 2 {\V hh \ 2 ) n -> as 
71 — > 00. That is, a non-decaying particle can be entirely 
converted into a light mass state and become unconfincd. 

For completeness of the discussion, we qualitatively 
consider two more cases. Their quantitative analysis is 
beyond the scope of this paper. In case (A), let's assume 
the box is completely sealed (no membranes), so both 
mass states are kept inside. Then an equilibrium shall 
exist. There are two regimes: coherent and incoherent. 
If the wave packets are wide enough (their widths depend 
on the properties of production/detection, reflection and 
scattering processes) , then they overlap and interfere co- 
herently. Scattering is also coherent in this regime and 
the transition matrix elements shall be oscillatory func- 
tions of space and/or time. So, we expect that the equi- 
librium state shall depend on masses, energies and mo- 
menta of mass states, the size of the box and, perhaps, 
parameters of the scattering potential as well. In the 
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second regime, when all wave packets are localized and 
their overlap is vanishing, the particle wave function is 
their incoherent sum. Then, scattering of mass states 
occurs independently, so an equilibrium is set by the de- 
tailed balance , \mh) T± \mi). Equality of the rates of 
the forward and reverse processes, anVih = aiVhi, and 
the normalization condition \a,i\ 2 = 1 uniquely deter- 
mine the equilibrium probabilities \di\ 2 , i.e., the mass and 
flavor composition of the system. In case (B), let's now 
put our original box with a particle (Fig. [T]) into a larger 
sealed box of size Lbox.2 > £box and consider the incoher- 
ent regime only (note that light states can come back into 
the smaller box, but heavy states cannot go out into the 
larger one). The mass conversion is now skewed toward 
\m h ) — > \mi) process, because the rate of the opposite 
process is reduced by a factor of Lbox/^box,2- Hence, 
the equilibrium composition is skewed towards the light 
states. In the extreme case of Lbox,2 ~> oo, our original 
result, la^ 00 '! 2 = 0, is recovered. 

So far we considered an idealized model. In a realis- 
tic physical model one shall solve the Dirac equation to 
obtain accurate amplitudes and transition probabilities 
in the four-dimensional space-time. Note that in some 
quantum setups the V and M processes cannot be taken 
as independent. The analysis of such cases is interesting 
but is subject to a dedicated study elsewhere. We also 
want to keep the discussion fairly general, hence consid- 
eration of fine retails or specific models is beyond the 
scope of the paper. Here we consider a case when scat- 
tering and propagation are effectively decoupled, so the 
V and M. processes are independent. 

The plane wave approximation is used for the V. A 
free particle solution of the Dirac equation for ith mass 
state propagating in a certain direction is ~ e lPiX , which 
can be normalized by the box volume; here P = (e, p) 
and X = (t, x) are four- vectors, and we use h = c = 
1 throughout. In the perturbation theory, the matrix 
elements are V 3l ~ (m 3 -| / d A X e~ lP 'i x V(x.)e tP * x |m 4 ), so 
that the time-dependent integral singles out into a 5- 
function, J dt e ie j*~ leit ~ S(e'j — e»), which is just the 
conservation of energy. Note that m; — > nih process can 
be kinematically suppressed if the primary (mi) is slow, 
i.e., when mf + pf < m 2 . In the m/j — > fn\ conversion, 
the secondary (mi) has a larger momentum, mf + p' 2 = 
m h + Ph ■ For a non-relativistic primary, the secondary 
can be relativistic and non-relativistic, respectively, 

p[ — m h and v[ — (2Am/m) 1/2 , (4) 

where ph was neglected, and in the latter case mh ~ mi = 
m and Am = mh — mi. 

We now consider the physics of the A4-process. What 
in a real world would play a role of a box with mem- 
branes? We suggest to consider a gravitational potential 
of any origin (for instance, a star or a black hole, a galaxy 
with its dark matter halo, a cluster of galaxies, etc.). 



Obviously, one can always adjust the masses and mo- 
menta of \mi) and |m^), so that the heavy state is bound 
(trapped) in the potential (vh is smaller than the escape 
velocity) and the light state is unbound (e.g., relativis- 
tic or non-relativistic but just fast enough). Since the 
gravitation potential is macroscopic, the evolution can 
be treated in the wave-packet approximation ^ for each 
mass state. For simplicity, we consider a non-relativistic 
one-dimensional packet, 

ip(x,t) = exp {ir] t (x - x t ) 2 + ipt(x - x t ) + i(j>t} , (5) 

where m, 4>t are complex, p t , Xt are real and all are 
functions of time. If (ip\ip) = 1, then, obviously, Xt — 
(tp\x\^) and p t — (ip\p\ip) take their simple classical 
meaning. In turn, r\ t determines the width of the wave 
packet and the phase <p t can be shown to be the action 
along the path. For a locally quadratic potential, U (x) = 
U + U'(x - x t ) + U"(x - x t ) 2 , i.e., with U'" and higher 
terms being identically zero, where Uq — U(x t ). l U' = 
d x U\ x=Xt , etc., it can be shown by direct substitution into 
the Schrodinger equation iip = Hip with the Hamiltonian 
H = -(2m)- 1 d 2 x + U(x) that 

x t = dH/dpt, p t = -dH/dxt, (6) 
fit = -2rff/m - U"/2, § t = ir\ t jm + x t p t - E, (7) 

where E = H(x t ,pt) = p 2 /2m + U(xt). Particularly 
important for us is the first pair, Eqs. ([6]), which demon- 
strates classical motion of the packet. It is fair to say 
that mass states propagate along geodesies. Whether a 
mass state is trapped or not is determined by it's velocity 
and the escape velocity of the potential. 

Implications — Below we offer a few possible scenar- 
ios where mass conversions can occur. A more detailed 
analysis is beyond the scope of the paper. 

The cosmic neutrino background (CNB) is the relic 
neutrinos [6] with the present day temperature of < 
2 K. Although neutrino masses are unknown, one can es- 
timate them from oscillations [7] to be m ~ (Am 2 ) 1 / 2 ~ 
0.05 — 0.007 eV, implying that cosmological neutrinos are 
nonrelativistic and the heaviest ones can be trapped in 
dark matter halos with escape velocities of order a few 
thousand km/s. Hence, if one ever be able to detect CNB, 
one can notice that it's composition is distorted from 
the the expected equipartition between the flavors. As- 
suming the non-relativistic cosmological neutrino cross- 
section cr„ ~ 10 -60 cm 2 (it can be greatly enhanced by 
coherent effects, however [6]), the mean density of nu- 
cleons in the halos n ~ 10 -3 cm" 3 , the thermal velocity 
of v ~ 10~ 2 , we estimate the number of scatterings to 
be ~ na^vtn ~ 10~ 47 per particle in the Hubble time, 
tn = 4 x 10 17 s. The effect is very tiny and unobservable 
(unless nonrelativistic neutrinos interact more strongly 
then we currently think) , though it may have some effect 
for the ultimate fate of the Universe, at tu oo. 

Another interesting possibility is axion-photon mix- 
ing, if axions constitute cold dark matter (CDM) jSJ[H]. 
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In this scenario, axions can scatter off cosmic mag- 
netic fields. Using an experimental upper limit on 
axion-photon coupling g ai ~ 10~ n GeV -1 , we esti- 
mate the interaction/conversion probability [5] to be ~ 
0.1( ffQ7 /10- u GeV _1 )(.B/l nG)(L/10 kpc). For galac- 
tic halos with the typical galactic field strength of, say, 
B ~ 1 — 3 fiG and the halo sizes L ~ 30 — 100 kpc, 
the probability is of order unity per passage. In galaxy 
clusters with B ~ 0.1 /iG and sizes of a couple of Mpc, 
the rate can be similar. Interestingly, this is the most 
optimal regime of conversions (one per passage) and the 
cosmological effect can be noticeable. If occurs, it can 
lead to substantial "evaporation" of axion CDM halos 
on the Hubble time-scale. If the effect is not observed, it 
can further constrain the coupling constant, g a ^. 

The third scenario deals with the CDM in the weakly 
interacting massive particle (WIMP) sector, which is still 
an attractive model [9] [10] . Usually, it involves a sta- 
ble lightest (tens of GeV-TeV-scale) supersymmetric par- 
ticle - supposedly a mixture of several supersymmetric 
particles, in fact. For high mass degeneracy, which we 
consider now, decay channels can be kinematically for- 
bidden, so more than one particle can be stable. Taking 
rather crudely Am ~ MeV and m ~ TeV, we estimate 
the velocity of a light state after mn — >• mi conversion, 
Eq. to be v[ ~ 300 km s _1 , which is a cosmolog- 

ically interesting number. It is smaller or comparable 
to the escape velocity from halos of large galaxies and 
galaxy clusters, but it is larger then v esc for halos of 
dwarf galaxies, which can be as low as a few tens km/s. 
A halo dichotomy is then expected: the number of small 
halos must be smaller than what CDM simulations pre- 
dict whereas the population of large halos is mainly unaf- 
fected. Such an effect can be an alternative explanation 
to the 'missing satellite problem' and, perhaps, it can 
also affect halo cores [TTJ. For the effect of conversions 
to be significant, the rate shall be at least a few reac- 
tions per Hubble time, which seem to rule out scattering 
off normal matter. Dark matter self-interaction can ac- 
commodate the needed large cross section, though sim- 
ple models seem to be disfavored [TT]. The models with 
the velocity-dependent cross-section and the Sommerfeld 
enhancement [5], which is the most profound at low ve- 
locities (that is, in dwarfs, again), shall be tested before 
this scenario is ruled out entirely. The disadvantage of 
this scenario is also in that it required strong degener- 
acy, which needs an explanation. However, a similar but 
ad hoc model with GeV-scale WIMPs with a keV-scale 



mass-degeneracy seems to explain some direct detection 
experimental data [12 . 



The author thanks Nima Arkani-Hamed and Peter 
Goldreich for thoughtful discussions, Masataka Fukugita 
for interesting suggestions and Freeman Dyson for en- 
couragement. The author acknowledges support from 
The Ambrose Monell Foundation (IAS) and The lb 
Henriksen Foundation (NBIA), as well as grants AST- 
0708213 (NSF), NNX-08AL39G (NASA) and DE-FG02- 
07ER54940 (DOE). 



edu 



[9 



[10 



mvm'iSias 

Also at: Niels Bohr International Academy, Blegdamsvej 
17, DK-2199 K0benhavn 0, Denmark; Department of 
Physics and Astronomy, University of Kansas, Lawrence, 
KS 66045; Institute for Nuclear Fusion, RRC "Kurchatov 
Institute", Moscow 123182, Russia 

B. Pontecorvo, Soviet Journal of Experimental and The- 
oretical Physics, 6, 429 (1958). 

S. Nussinov, Physics Letters B, 63, 201 (1976); 

B. Kayser, Phys. Rev. D, 24, 110 (1981). 

C. Giunti, C. W. Kim, and U. W. Lee, Phys. Rev. D, 
45, 2414 (1992); C. Giunti, C. W. Kim, J. A. Lee, and 
U. W. Lee, ibid., 48, 4310 ( 1993) , |arXiV :hep-ph/9305276l 
J. T. Peltoniemi and V. Sipilainen, Journal of High 
Energy Physics, 6, 11 (2000 ) |arXiv:hep-p h/0004162 
C. Giunti, ibid., 11, 17 (2002), |arXTv7he p-ph/0205014, 
E. W. Kolb and M. S. Turner, The early universe, Fron- 
tiers in Physics (Reading, MA: Addison- Wesley, 1990). 
A. D. Dolgov,|Physics of Atomic Nuclei, 71, 2152 (2008)| 
larXiv:0803. 38871 

S. M. Bilenky, C. Giunti, J. A Grifols, and E. Mass o, 
|Phys. Reports, 379, 69 (2003)| |arXiv:hep-ph/0211462| 
J. E. Kim and G. Carosi, Reviews of Modern Physics, 82, 
557 (2010)] N. Bassa n A. Mirizzi, an d M. Roncadelli, 
ArXiv e-prints (2010), [arXiv:1001.5267| 
N. Arkani-Hame d, D. P. Finkbeiner, T. R. S latyer, 
and N. Weiner |Phys. Rev. D, 79, 015014 (2009)] 

larXiv:0810.07l3l 

J. L. Feng, Journal of Physics G Nuclear Physics, 32, 1 



[12] 



(2006), arXiv:astro-ph/051104 3) 

~K. Kravtsov, |Advances in Astronomy, 2010 (2010) 



doi:10.1155/2010/281913, arXiv:0906.3295; R. Kuzio 
de Naray, G. D. Martinez, J. S. Bullock, and 

J. Lett 



M. Kaplinghat, Astrophys 
larXiv:0912. 35181 
R. Essig, J. Kaplan, P. Schuster 
e-prints (2010), |arXiv:1004.069l| 



., 710, L161 ( 2010) 
and N. Toro, ArXiv 



